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Abstract 

Approximation of the previously estimated excess specific heat (J^^'^'^^^ /T^ 
of two fluoropolymers, PTFE and PCTFE, is presented using Orthonormal 
Polynomial Expansion Method (OPEM). The new type of weighting functions 
in OPEM involves the experimental errors in every point of the studied ther- 
mal characteristic. The investigated temperature dependence of the function 
(jexcess ^rp5 jg described in the whole temperature ranges 0.4-^8A' and 2.5-^7K 
respectively for PTFE and PCTFE as well as in two subintervals (0.4 2)K, 
(2.5 -7- 8)K for PTFE. Numerical results of the deviations between the eval- 
uated (7e2;cess^ji5 ^g^ig^ g^j^d their approximating values are given. The usual 
polynomial coefficients obtained by orthonormal ones in our OPEM approach 
and the calculated in every point absolute, relative and specific sensitivities of 
the studied thermal characteristic are proposed too. The approximation pa- 
rameters of this type thermal characteristic are shown in Figures and Tables. 

Key words: orthonormal and usual polynomial approximation, low-temperature ex- 
cess specific heat of two fiuoropolymers - PTFE and PCTFE 

I.INTRODUCTION 

The unusual thermal properties concerning the low-temperature specific features 
of the heat capacities of the polytetrafiuoroethylene (PTFE) and polyclorothrifiuo- 
roethylene (PCTFE) were considered in an earlier paper [1]. The estimated there 
excess specific heat over Debye contribution below 10 K of these fiuoroplasts was 
discussed in the frame of the recently developed Soft-Potential Model (SPM). 
The present study is devoted to the mathematical description of the low-temperature 



excess specific heat of these semi-crystaUine polymers applying Orthonormal Poly- 
nomial Expansion Method (OPEM) [2]. 

lI.EXCESS SPECIFIC HEAT DATA 

The previous work [1] clarifies several points concerning the specific heat pe- 
culiarities observed in Cp/T^ temperature dependencies of PTFE and PCTFE, a 
maximum appearing around T^^x— 4 K for both polymers and a found shallow 
minimum for PTFE centered at Tmm — ^ K. Following the Soft-Potential Model 
(SPM) [3, 4, 5, 6], supposing a coexistence of acoustic phonons with quasi-locahzed 
low-frequency (soft) modes in glasses and successfully applied by us to some chalco- 
genide glasses [7, 8], the specific heat data, taken from the available calorimetric 
measurements at low-temperatures (T < 10 K) for the studied fluoroplasts [9, 10], 
were described in a paper [1] . The Cp components are: 

i) Cj^*^ - a linear contribution, described by double-well potentials, conditioned by 

the thermal exitations of the tunneling state (TLS); for PTFE it was established to 
predominate at T < 0.2 K. 

ii) (C'^)«cou _ cubic Debye contribution, evaluated and discussed in details in a 
paper [1]. 

iii) C'ea^cCSM) _ excess specific heat (a soft mode) contribution of the quasi-harmonic 

exitations, described by single-well potentials. 

(jexc{SM) component was evaluated [1] by difference between the measured specific 
heat Cp [9, 10] and the sum (Cj^^ + (C^)"^) as follows. 

The temperature dependencies of the three specific heat components in Eq. 1 are: 

The abbreviations of Two Level State, Debye and excess specific heat (Soft Modes) 
are marked with TLS, D, and exc{SM), respectively. Here, Ctls was determined 
by the Cp experimental data of Nittke et al [9]; (Cd)""^°", the true elastic coefficient, 
was calculated [1] by the macroscopic parameters of the investigated materials and 
the average sound velocity Vs-, evaluated [1] for PTFE (PCTFE) by available mea- 
surements of transversely and longitudinally polarized 10 (5) MHz ordinary sound 



both polymers changes its value about 0.3 % within and 10 K. This fact allows 

us to accept up to 10 K a constant Debye coefficient of the Debye specific heat 
contribution by acoustic measurements . 

The quantities Cd up to 10 K and Ctls are independent of the temperature. In 
accordance with the SPM, the softening of the lattice vibrations leads to increasing 
of the density of states of the harmonic oscillators with rising of their energy E. This 
increase appears to be proportional to the energy as E'^ leading to proportional to 
term in Cp{T) at enough low-temperatures, i.e. the quantity (^6x0(51^) temper- 
ature independent only in a narrow low-temperature range, < T < 1.5T^j„. It is 
worth to mention that the quantity C^M^^) evaluated for some glasses [11, 12, 13] 
confirms this prediction of the SPM. For the studied polymers [1] T^m was esti- 
mated to be very close to 1 K. But in wider low-temperature range the temperature 
dependencies of the quantities Cp^^^^^^T^ and the excess specific heat C'pa;c(S'M) 
presented, respectively in Fig.l and in the inset of the Fig.l, are temperature de- 
pendent. As it can be seen in the inset of Fig.l, the Soft Modes (SM) delocahzed 
at T = 8 K for PTFE and at T = 7 K for PCTFE. It is important to note that 
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Figure 1: Temperature dependencies of the excess specific heat presented as C^^'^/T^; 
inset: Temperature dependencies of the excess specific heat C^^"^ [mJ/gK] of the 
PTFE and PCTFE 



(jexcjrj^b function changes its temperature behavior just at T = {2>/2)Tmin- 

III. MATHEMATICAL APPROACH 

Our Orthonormal Polynomial Expansion Method(OPEM) and its applications 
in cryogenic thermometry are presented in papers [2, 14, 15, 16]. Some important 
features of OPEM concerning cryogenic thermometry at the approximation of ther- 
mometric characteristics of different type low- temperature sensors are protected by 
a patent for an invention [15]. Our OPEM is defined on Forsythe [16] three-term re- 
lation for constructing orthogonal polynomials over discrete point set with arbitrary 
weights in the term of the least square method. The one-dimensional recurrence 
for generation of orthonormal polynomials A; = 1, 2, . . .} and their derivatives 

{^P,m = 0, 1, 2, . . .}, in OPEM is: 

= 7.+i[(? - - (1 - 5.o)/3.*K(g) + m^>t-'\q)] (3) 

OPEM is a development of the Forsythe approach for receiving derivatives and 
integrals with fourth term in the Eq. (3). The polynomials {^1°''} satisfy the 
orthogonality relations: 

M 

i=l 

over the point set {qi,i — 1,2, .. . M} with weights Wi — l/af , depending on errors 
(Tj in every point. The approximating values of the function and its m-th 

derivative /M'^PP'",{m — 0, 1, ..} are calculated by 

yMapp.(^) ^ ^ a,^^;^\q) = cm'- (4) 

k=0 k=0 

The optimal degree N of the approximating polynomials in Eq.(4) is selected by the 
algorithm, combining the following two criteria. First, the fitting curve should lie in 
the error corridor of the dependent variable {qj, /J^^ ±aj,j= 1, ...M). 

Second, the minimum should be reached. 

M 

E ^^Afr - /D -N-1)^ ruin. (6) 
When the first criterion is satisfied, the search of the minimum stops. The 



involving the errors in independent variable was published in a paper [17]. The last 
version with obtaining of usual Ck coefficients from orthogonal ones au from Eq. 4 is 
developed in our work, RSI-2005 [2]. 

IV. APPROXIMATION RESULTS 
A. Orthonormal expansion 

The temperature dependence of the C^^'^jT^ function is described by orthonor- 
mal polynomials in the whole temperature range 0.4 8 \K\ or in two subinter- 
vals 0.4 -^ 2 \K\i 2.5 -^ 8 \K\ for PTFE, and in the temperature range 2.5 -^ 
7 \K\ for PCTFE, using the new type of weights, Vl^c-'^Vr^. The studied ther- 
mal characteristic C^^'^/T^ and its sensitivities, absolute d{C^^'^/T^)/dT, relative 
[l/(C^^7T5)]d(C^^7T5)/dT and specific, [T/{C'''"'/T^)]d{C^'"'/T^)/dT, evaluated 
in every point, are shown in linear-linear plot in Figs. 2 and 3 for PTFE and PCTFE, 
respectively. It is worth to note that the subintervals for PTFE are chosen by the 
temperature behavior of the specific sensitivity of this polymer (see Fig.2). 



8 



o 
c 

3 



u 
o 

(A 
0) 



(A 
C 
0) 



4- 
2 

o^ 

-2 
-4 
-6 
-8 
-10 



O (1/C^'"')dC^'"'/dT 
^ (T/C*''°)dC*'"'/dT 




012345678 

Temperature 



PTFE 



o ooooo 



2 3 4 5 6 

Temperature 



o 



4 

2 



-2 

-4 

-6 

-8 

-10 



8 



Figure 2: Temperature dependencies of the relative [l/{C'^^'^/T^)]d{C'^^'^/T^)/dT 
[K-'^] and specific [T/{C^'"'/T^)\d{C^'^''/T^)/dT [-] sensitivities of the PTFE; inset: 
Temperature dependencies of the absolute sensitivity d{C^^'^/T^)/dT. 
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Figure 3: Temperature dependencies of the relative [l/(C^''yT^)]d(C^''yT^)/dT 
[K-^] and specific [T /{C^'^^ /T^)]d{C^''^ /T^)/ dT [-] sensitivities of the PCTFE; in- 
set: Temperature dependencies of the absolute sensitivity d{C^^^/T^)/dT [mJ/gK'^]. 

By definition the weighting function W'-''""'/'^^ is l/cr^, where is a variance 
of the thermal characteristic C'^^'^/T^^ versus temperature T. In our investigation 
this variance is accepted to be, correspondingly square of the absolute heat capacity 
resolution /^C^, determined by the experimental specific heat accuracy as follows: 
(AC^^;^)i = 0.025(C^^7r^)i [mJ/gK^\ for the first approximating interval of PTFE 
and (AQ^^^)i = 0.1(C^^7T5)i [mJ/gK^] for both, the whole PCTFE approximating 
interval and the second approximating interval of PTFE. Here the weights, W^"'"'^'^^ 
are expressed by the relations: 

^^,c-7T^). = 1/(AC^-/T„5,J, = 1600/(C^^7T5)2 [mJ/gK^]-^ (7) 

for the first approximating interval of PTFE, or 

(M^^-/^^), = l/{^Cli:), = 100/(C-7T^),^ [mJ/gK^^ (8) 

for the PCTFE and the second approximating interval of PTFE. The deviations 
A(C^^'^/T^)j between experimental and approximating values of the excess specific 
heat are estimated in every point by the expression: 



The temperature behavior of the calculated differences A{C^^'^/T^)i, the root mean 
square deviations RMS'" and the mean absolute deviations MAD'~^ respectively for 
PTFE and PCTFE excess specific heat approximations are shown in the Figs. 4, 5 
and 6. The RMS deviation is more popular. It is given in another our paper [2]. 
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Figure 4: Temperature dependencies of A(C^^7r^), (A(C^^7r5)^^^^ MAD^"'/^' 
and RM S^"'"' /'^^ of the PTFE for the whole approximated temperature 
range;inset:Temperature dependence of the weighting function w{T) for FTPE . 

The characteristics MAD is defined as follows: 

MAD^ ^^Yl I A(c^^7r^)i - A(c«^vr5)| , (lo) 

i 

where A{C^^^/T^ = ig A(C^^7T5)i and A(C^^7T5)i are calculated from Eq.(9). 
Following the cited criteria in Eqs.(5,6) the deviations A{C'^^'^/T^)i are in the error 
corridor (see Eq.5). As a result of our approximation, the polynomial degree for 
temperature range of PTFE subintervals is lower than for the whole temperature 
range. For C^^^/T^ vs. T approximation the optimal degree N as well as some main 
characteristics including the overall approximation characteristics: RMS and MAD 
of AC^^'^/T^ vs. T approximation, the weighting function W'-'' and a goodness of fit 
are presented in the TABLE 1. 

B. Usual expansion obtained by orhonormal one 

The next step in our approximation is done here. The OPEM is extended by 
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Figure 6: Temperature dependencies of A{C^^yT^), A(C^^'^)„cr, MAL)^'"'/^" and 
RM S^""^" of the PCTFE for the whole approximated temperature range. 



Table 1: QPEM approximations of C^^^<^'' jT^ for PTFE and PCTFE 



T 




C 


N 






RMS^ 




MAD^ 




x' 


[K] 




[liJ/gK^] 






[mJ/gK^] 


[mJ/gK^ 










0.4^ 


2 


152^1 


6 


6.91 ; 


< ^ 1.60 X 10^ 


0.121 X 10" 


-2 


0.828 X 10" 


-3 


1.38 


2.5^ 


8 


1^0.015 


4 


2.65 > 


; 10® ^ 4.44 X 10^^ 


0.265 X 10^ 


-4 


0.208 X 10" 


-4 


0.57 


0.4^ 


8 


152^0.015 


8 


6.91 X 10"^ ^ 4.44 X 10^^ 


0.109 X 10" 


-2 


0.326 X 10" 


-3 


0.95 


2.5^ 


7 


1.75 0.031 


6 


3.26 X 10'^ 1.04 X 10^^ 


0.428 X 10" 


-4 


0.306 X 10- 


-4 


0.58 



Table 2: Usual coefficients for (^e^cess^j^s j^j^Q ^faJ/gK^] vs. T approximation 



material 


PTFE 


PCTFE 


Raiuj(T{K) 


0. 1 ^ 2 


2.5 ^ 8 


0. 1 ^ 8 


2.5 ^ 7 


RangeiC^^^/T^) 


152 1 


l-^ 0.015 


152 -=-0.015 


1.75 -^ 0.031 


Co 


.022854 


.003118 


.016067 


.182459 


Cl 


-.001420 


-.009158 


-.063937 


-.667996 


C2 


.105876 


.011825 


.103755 


.997805 


C3 


-.582978 


-.006389 


-.070355 


-.769121 


C4 


-.390617 


.001190 


-.008661 


.320591 


C5 


.845833 




.048421 


-.068068 


C6 


.696453 




-.033867 


.005687 


C7 






.010411 




C8 






-.001230 





is carried out for the investigated intervals. The calculations are made for the 
(jexcjrpb _ f[T^ descriptions, in two runs: first - in an interval [-1,1], and second - 
in the input intervals. Optimal values of usual polynomial degrees are chosen using 
two new criteria. 

The first criterion is: 

max mm, (11) 

where f"-PP^'^ is the approximating function defined with usual expansion. The sec- 
ond criterion is: 

A{cfe}f ^ min, (12) 

where {Acfe} are inherited errors in usual coefficients, defined in OPEM, discussed in 
a paper [18]. In the TABLE 2 the usual coefficients for C^^'^/T^ vs. T approximation 
are presented. 



In conclusion, the heat capacity components of the studied fluoroplasts are well 
defined earlier [f] using the SPM. In this study, the temperature dependencies of 
the excess specific heat of PTFE and PCTFE below 10 K presented as C'^'^^^^T^ 
are described by an Orthonormal Polynomial Expansion Method (OPEM). An ap- 
proximation of the above mentioned thermal characteristic by usual polynomials, 
obtained by orthonormal ones, is given too. The numerical results for approxima- 
tion characteristics as root mean square deviation, mean absolute deviation and 
normalized x square assure good accuracy for calculated values of physical quantity 
in the input temperature intervals. Our approach for description of this characteris- 
tic proposes possibility for the future analyzing the low-temperature excess specific 
heat of the materials with glass-like behavior. 
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